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Abstract 

We consider the effects of higherloop corrections to a Schwinger-Dyson 
equations for propagators. This is made possible by the efficiency of the 
■ methods we developed in preceding works, still using the supersymmetric 

Wess-Zumino model as a laboratory. We obtain the dominant contribu- 
tions of the three and four loop primitive divergences at high order in 
\ perturbation theory, without the need for their full evaluations. Our main 

£S| ■ conclusion is that the asymptotic behavior of the perturbative series of the 

renormalization function remains unchanged, and we conjecture that this 
will remain the case for all finite order corrections. 
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1 Introduction 



In preceding works [U[2j[3] , we have explored the solution of a simple Schwinger- 
Dyson equations for propagators, taking as a laboratory the supersymmetric 
^ \ Wess-Zumino model. In this way, the perturbative series of the renormaliza- 

tion group functions could be computed up to the high order of 200 loops and its 
asymptotic behavior predicted. However, these results have been obtained for 
a Schwinger-Dyson equation which includes only the one loop primitive diver- 
gence: the full Schwinger-Dyson equation includes an infinite number of higher 
loop primitive divergences. We here explore the effect of the first additional 
terms of the Schwinger-Dyson equation in the supersymmetric Wess-Zumino 
model: the three-loop and the four-loop terms (there are no primitive two-loop 
divergences in this model). This is by no means a simple task, since we cannot 
be content with the computation of the diagrams, but need to obtain a Mellin 
transform with eight variables in the three-loop case and eleven in the four-loop 
case. Furthermore, even if there is only one diagram at the level of supergraphs 
for each of these orders, they correspond to a sum of differing terms with pos- 
sibly fermions. A scalar theory would present its own difficulties, since in this 
case, there will be vertex corrections and overlapping divergences which present 
conceptual problems. 
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We can however obtain valuable answers, through a chain of simple steps. 
The evaluation of the Mellin transform can be made with the scale invariant 
vacuum diagrams obtained by connecting the two external legs of the diagrams. 
The high degree of symmetry of these diagrams, the ^£3,3 complete bipartite 
diagram at three loops and the cube at four loops, simplify many computations. 
We concentrate on terms which have the potential to change the asymptotic 
behavior of the perturbative series, which is not the case for single terms in the 
Taylor expansion of the Mellin transform. The dominant contributions there- 
fore stem from the leading parts of the asymptotic expansion of the Taylor 
series, which can be deduced from the singularities near the origin of the Mellin 
transform. These singularities in turn stem from the divergences which appear 
when some subdiagrams become scale invariant. The poles are easy to deter- 
mine and the values of the residues which dominate our computation reduce to 
known bosonic integrals. 

We shall first recall some basic points of our precedent works, which motivate 
the present study. We next present the link between the potential divergences of 
subdiagrams and the poles of the Mellin transform, allowing a complete study 
of the associated residues. For simplicity, the presentation is made for scalar 
diagrams, but the possible numerators, coming from fermion contributions in 
our case, add only inessential computational complexities. We then consider 
the numerators for the combination of diagrams imposed by supersymmetry. 
Their full evaluation in a parametric representation of the diagrams is beyond 
the possibilities of the classical methods, if we want to keep the high degree of 
symmetry which one expects: we therefore limit ourselves to the cases which 
are necessary for the evaluation of the dominant residues. We then combine 
all these results to obtain the leading asymptotic contributions of these two 
diagrams. 

2 Summary of previous work. 

Our solution strategy for the Schwinger-Dyson equation was developed in for- 
mer works [U [2] , based on the proposal of [I] . A first ingredient is a renormal- 
ization group argument which allows to obtain a full propagator as a series in 
the logarithm of the norm of the momentum from the renormalization group 
functions for a minimal computational cost. The power of methods using Hopf 
algebras of graphs [5j EJ [7] allows to show that these renormalization group 
arguments apply to the solutions of truncated Schwinger-Dyson equations. It 
has further been remarked in j3j that not only the different derivatives of the 
propagator can be computed, but also some infinite sums of derivatives. 

In this simple case with only renormalized propagators, the remaining prob- 
lem is the evaluation of primitively divergent diagrams with full perturbative 
propagators. Since perturbative propagators are polynomials in the logarithm 
of the squared momentum, all such evaluations can be derived from a sin- 
gle analytic function, with parameters appearing as exponents of the different 
propagators of the diagram. Taylor expansion then allows to obtain the effect 
of any power of the logarithm of the propagators. The use of such parameters is 
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the basis of the analytic regularization which enjoyed a brief popularity before 
being superseded by dimensional regularization [8j [U [10] . Such Mellin trans- 
forms have now a rich history as a tool for the evaluation of Feynman diagrams: 
in particular, Mellin transforms are instrumental in [11] to compute multiloop 
diagrams from single loop parts. 

What will be of tremendous use in what follows is that the Mellin transform, 
which is defined a priori as an analytic function on a domain of convergence 
depending only on the real parts of the parameters, can be extended on the 
whole complex domain as a meromorphic function. Its poles are linked to di- 
vergences, when subdiagrams become scale invariant. In a parametric integral 
representation of these Mellin transforms, factors of the inverse T function ap- 
pear, which give a number of zeros. For even integer dimensions and integer 
values of a number of the parameters, these zeros cancel but a finite number of 
poles, so that one obtains a rational function. Such simplifications will allow to 
give dominant contribution from a priori very complicated diagrams. 

3 Scale invariant diagrams. 

We begin by the study of bosonic diagrams, avoiding the complexities stemming 
from numerators in the super symmetric case. We start from a diagram with 
two external legs G, to be evaluated with an external momentum p, with a 
Mellin transform variable /3j associated to each edga3- 

It is convenient to use a parametric representation of the diagram G of loop 
number L, with Schwinger parameters a^, integrated from to oo, for each edge 
i. The first and second graph polynomials U and V, respectively of degrees L 
and L + 1, allow for a simple expression of the integral on loop momenta in 
dimension D: 

Usually one simply integrates on the overall scale of the a^, but we can introduce 
an additional variable ao to express the product of powers of U and V in terms 
of U = V + clqU, which is the polynomial associated to the vacuum diagram 
G obtained by completing the original diagram with an edge joining the two 
vertices linked to the exterior. With the exponent /?o such that the sum of 
the fa is (L + l)D/2 for a diagram with L loops, takes the following 

very symmetric form, depending only on G, apart from the dependence on the 
exterior momentum: 

W] = ^ D/2 jlYW) I Xiida^U-VH&ai ~ ^ (2) 

where now the products include factors with index i = associated to the edge 
closing the graph. The values of the give an overall scale invariance to the 

1 Physicist are used to call them lines, but we shall stick to the terminology of graph theory. 
Line will be used to denote a path formed by Fermionic propagators. In the supersymmetric 
cases, we have to sum on the different configuration of Fermionic lines. 
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integrand in this formula, which is fixed by the 5 factor. The sum is primed in 
this 5 factor to indicate that any subset of the dj can be included in this sum, 
without changing the value of the integral. This freedom makes the analysis of 
the singularities of this integral easier. 

Consider indeed any 1PI subdiagram K of G and the new diagram G/K 
defined through the contraction of the subdiagram K, i.e., deleting the edges of 
G belonging to K and identifying their ends. If one gives weights 1 to the edges 
of K and weights to those of G/K, the lowest weight terms of Uq factorize 
into Uq/kUk [12] • If we introduce a new variable a appearing as a factor in all 
the variables associated to K, the integral in equation ([2]) becomes, if K has I 
loops, 

/ ^(da^^diZ'^-mZ'a^-l) J daa^ 1 (a"U G/K U K +- ■ -)- D / 2 , (3) 

with (3 = Y^i'eK@i'- The indexes i' are associated to the edges in K and the 
i" to the edges in G/K. The second 5 is there to compensate the introduction 
of the additional integration variable a. The a integral is convergent in when 
f3—£D/2 is positive and can be extended as a meromorphic function with simple 
poles when /? — ID/2 is or a negative integer, with residues proportional to 
the integrand and its derivatives at a = 0. For f3 — £D/2 = 0, the integrations 
over the a» associated to K and G/K factorize, so that the residue of is 
given by 

^ p _ w/2=0 l[(3 i ] = l/T{D/2)l K [^\l G/K [p i ,,], (4) 

The condition on /3 is exactly what is necessary for both G/K and K to be 
scale invariant, so that the integrals are well defined. For the further poles 
(3—£D/2 = —k, the denominator is still factorized, of the form U K D l U G ^ 2 1 , 
with I an integer less or equal to k. The residue therefore is a sum of factorized 
integrals. 

We conjecture that these poles are the only singularities of the meromorphic 
extension of with double and higher poles associated to chains of inclu- 

sions of diagrams G\ C G 2 C • • • C Gi = G. The combinatorics associated to 
the different decompositions in chain of subdiagrams should be easy to control 
with the core Hopf algebra introduced in \13\ I14j. This algebra is similar to the 
Connes-Kreimer Hopf algebras of diagrams [6], but without the restriction to 
diverging subdiagrams in the left part of the coproduct. The quotient diagram 
in the right part of the coproduct will therefore have vertices with any valence. 
Indeed, with a proper choice of the variables any subdiagram can become 
divergent. The order n poles of should be read from A n 

gE For a L loop 

diagram, the residues of the poles of order L — 1 come from chains of inclu- 
sions such that each quotient Gi + \/Gi is a one loop diagram: the corresponding 
integrals will be of the form 

J n(^KxE a *r D/2 ^K - 1), (5) 

2 A is used to denote the reduced coproduct, which is obtained by eliminating the trivial 
1 ® G and G ® 1 terms in the coproduct. 
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with P a polynomial whose degree is completely fixed by the requirement of 
scale invariance. For each monomial of P, this generalized B function evaluates 
to a product of Euler V functions of arguments shifted by integers, divided 
by r(— D/2 — I). Divided by the product of r(/3j) in the definition of I, this 
simplifies to polynomials in the /3j, with rational coefficients when D/2 is an 
integer. If there are no numerators, as is the case for the poles at the limit of 
the convergence domain, the residue is constant. 

4 Supersymmetry and numerators. 

As advanced in the introduction, we intend to consider a supersymmetric Wess- 
Zumino model defined in terms of a single chiral superfield composed of a com- 
plex scalar, a chiral fermion and a complex auxiliary field. At the order of 
our computations, the only consequence of the chirality of the fermion field is 
to divide by 2 the trace over spinor indices involved in fermion loops. Super- 
symmetry introduces both simplifications and complexities. The presence of 
fermions introduces numerators in the Feynman integrals. The special com- 
binations suppress a number of divergences, making the whole subdivergence 
structure simpler. On the other hand, these numerators which result from the 
combinations of a large number of terms have no simple general expression. In 
the cases of the three- and four-loops corrections however, the simplest diver- 
gences are given by purely bosonic diagrams. 

The only primitive three-loop correction in the Wess-Zumino model is given 
by the following non-planar diagram: 
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Indeed, for zero external momentum in the auxiliary field propagator, the 
combination of numerators becomes simply the product of the squared mo- 
menta on edges 1 and 7, so that the residue of the diagram is given by the 
tetrahedron diagram with exponents (/3i + jii — 1, /?3, /?4, f3e, fi-j + /3§ — 1). 
The same result can be obtained for the fermionic and bosonic propagators, 
but the necessary algebra is better done with computer assistance: we used 
Form from J. A.M. Vermaseren |15j . 

The first infrared divergences for any internal edge are likewise given by 
purely bosonic diagrams. The first infrared divergence indeed comes from an 
edge with a boson propagator, since the numerator in either fermionic or auxil- 
iary propagators would lower the infrared divergence degree. This can entirely 
fix the particle assignments on all edges, when we consider the edge 3 in the 
propagator for the auxiliary field, or strongly constrain the possibilities. In 
any case, the two edges meeting a zero momentum edge of bosonic charac- 
ter will share the same momentum q: either one of them is an auxiliary field 
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which contribute a numerator q 2 or the two are fermionic lines with numerator 
4 4 = Q 2 - In any cases, the two lines combine to a single bosonic line with 
exponent /3 + ft' — 1. 

For the divergence associated to Ps = 2, we simply obtain in the auxiliary 
field case the two- loop diagram with exponents (Pi + /3 5 — 1, p 2 , P4, Ps, P% + 
Pi — 1). This can be completed to a tetrahedral scale invariant diagram with 
an edge of degree Pq such that the sum of the Pi is 10 = AD/ 2 + 2. For the 
case of the fermion propagator, only three of the eight possible configurations 
survive when we impose that edge 3 is a boson propagator and they yield 
respectively factors of $2, —$4,$8 apart from the p\p\ factor common to all (all 
momenta are supposed to go from left to right). These factors add up to / 
p, so that supersymmetry is verified and the divergence is given by the same 
bosonic diagram as for the auxiliary field. In the bosonic cases, there are six 
contributions: the diagrams with only auxiliary and bosonic fields yield factors 
p 2 , p\ and p\ while the fermionic loops yield factors of 2p 2 ■ ps, —2p2'Pi and 
— 2p4-ps which add up top 2 , again trivially verifying supersymmetry. Combined 
with the symmetries of the diagram, this solves the cases of all the edges which 
do not meet the exterior ones. 

In the case of the divergence on the edge 1 for the propagator of the auxiliary 
field, the three possible field configurations contribute and they give factors 
p\p% p\p\ and 2p\ p^ ■ p§ for the fermionic loop. They add up to p\p 2 '. In 
the cases of fermionic and bosonic terms, it is easy to see that one obtains the 
same numerator, multiplied either by #f or p 2 . The strong symmetry between 
the different positions of the divergence suggests that, as in the purely bosonic 
highly symmetric formulation is possible, based on the complete bipartite 
graph -^3,3 obtained by closing the graph. In every cases, the divergence is given 
by the residue of the wheel with three spokes, which is known to be proportional 
to C(3). 

We now turn to the more challenging case of the four loop graph 
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The simplest situation is when we consider a divergence on the edge 5 for the 
auxiliary field: then the propagators 3 and 7 must be auxiliary fields and the 
numerator is completely fixed with the auxiliary field in 6. Apart from the 
factors coming from the neighbors of 5, we have the single factor p\. The cases 
of fermionic and bosonic propagators are more complex but can be easily solved. 
In the fermionic case, the fermion line can go either through 11 or 6. In the 
former case, there is only one possibility, which gives the numerator p\i>i\- In 
the latter case, one obtains (jfe +^4)^6^8 +#10), which gives ^6#fe^6 = p\i%- 
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The two add to form the result predicted by supersymmetry, p^p 1 . Likewise, 
one can show that for the bosonic propagator, one obtains a term p\ x p\ when 
11 is an auxiliary field, 2pu ■ p%p\ when 11 is fermionic and p\p\ in the last 
case where 11 is bosonic, for a total which is p\p 2 ■ 

With the infrared divergence on the edge 11, one obtains p\+2p^-p^+p^ = p 2 
with the auxiliary field, ($i +^2)^5 +$ta)(^9 +^10) = i>i>i> hi the fermionic case 
and [p\ + P2 + 2pi • P2){p\ + p\q + 2pg • p\o) = (p 2 ) 2 • Again, supersymmetry is 
simply verified. One remarks that the overall divergence degree is lowered. For 
a divergence on edge 3, one obtains a numerator proportional to p 2 and again 
supersymmetry relations. 

Let us detail the remaining case of a divergence on edge 1, which subsumes 
all necessary procedures. In the case of the auxiliary field, we must separately 
consider the cases where the propagator 6 correspond to the three possible 
field types. With a boson propagator, all other propagators are fixed and we 
get p\p\ beyond the common p 2 . With an auxiliary field on edge 6, three 
configurations are possible, which give terms p\p\\, p\p 2 and SPqPtPu, adding 
up to p1p\. Finally, if the propagator 6 is fermionic, the Fermion line can 
close either through edge 5 or 11, giving trfjfujfg^ej^) an d tr^^^i)- Now, 
P7 + P11 = P8, so that these two terms give 2p^ -p§p\. This three different terms 
have pi as a common factor, and the remaining factors add to p\ = p 2 . 

In the case of the fermionic propagator, the fermionic line goes first through 
edge 2 and we have four possibilities according to the next step being edge 4 or 6 
and the final one being edge 9 or 10. In each of these cases, a factor p\ appears, 
either directly when edge 8 is an auxiliary field propagator, or through the com- 
bination $$($7 +i>u) or through the combination of the fermionic loop (3, 5, 7, 11 
and the two possible auxiliary field positions on this loop. One therefore ob- 
tains ^§^2(^4 +$ta)(^9 +^10) = p\p 2 i>- I n the case of the bosonic propagator, a 
factor p\ = p 2 comes from the edge 2 which must be an auxiliary field and the 
three possibilities, auxiliary field on 9, auxiliary field on 10 and fermionic line 
going through both of them, have a p§ factor. The additional factors are p 2 ,, p 2 
and 2pg ■ pio which add up to p 2 . Therefore, also in this case is supersymmetry 
simply verified and the divergence is given by a simple bosonic diagram with 
exponents (ft - 1, ft + ft - 1, ft + ft) - 1, ft). 

Finally, we must consider the overall divergence of the diagram, which can 
be inferred in the limit of a vanishing exterior momentum. In the case of 
the auxiliary field propagator, it can be shown that the different numerators 
combine to give the product p\]p\p\, so that the residue is given by a scalar 
diagram with exponents (ft + ft — l,ft + fto — l,fti — 1), with all other 
exponents unchanged. The supersymmetry relations become harder to verify 
in this case, because they can no longer be verified at the level of the integrands 
through automated expansion. However, since the differences are proportional 
to the antisymmetric Levi-Civita tensor, the momentum integration is trivially 
zero. 

All these explicit computations tend to sustain the conjecture that the su- 
persymmetric Mellin transform can be obtained from a fully symmetric para- 
metric integration on the completed diagram, which is a cube. In any case, the 
singularities are always given by the wheel with four spokes graph, a diagram 
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whose evaluation is known to be proportional to £(5) (see for instance |16j). 

It should be noted that these higher loop primitive parts have contribu- 
tion of higher weight than the ones obtained through the computation of our 
preceding works. Indeed, in [I], the first term proportional to ("(3) appears at 
order four and the one proportional to £(5) at order six, at respectively one 
and two orders beyond the first contributions of the three and four loop primi- 
tives. This is consistent with the phenomenon of weight drop for diagrams with 
sub divergences described in [T7] . 



5 Applications. 

The cases we considered in the preceding section exhibit but a small part of the 
possible singularities, since they only are the first ones associated to subgraphs 
with one loop quotients. In addition to the further poles associated to the 
same subdiagrams, there are singularities associated to the subdiagrams with 
any number of loops and with any combination of these elementary singular- 
ities. However, these additional terms give subdominant contributions to the 
perturbative series of the (3— function and the limited knowledge we have will 
be sufficient for our purpose. 

The study of section [3] shows that the divergences are poles with denomi- 
nators built from the sum of the exponents associated to a subdiagram or its 
complement, since the exponent associated to the external propagator must be 
considered dependent variable. 

As in the fourth section of [3], where such general poles of the Mellin trans- 
form are studied, we consider the pairing as given by replacing the variables 
Xi such that = 1 — Xi by a derivative with respect to the logarithm of the 
momentum Lj. If we call Z(k,p) the object associated to the Mellin transform 



N(x,. 



k + Y7i=l X i 

it will be given by the following expression 

1 n 
l(k,p) = N(d Ll ,...,d Ln )l[G(L 

+ 2^j=l °Lj j=l 



(8) 



(9) 

-=L n =0 



The derivatives with respect to the Lj's can be converted through the use of 
the renormalization group to the expression 7 + j3 ad a . In the absence of vertex 
renormalization, j3 is simply 37 and the differential operator in a which is now 
in the denominator can be brought to the left hand side to give 

(k + j{p + 3ad a ))l(k,p)=N , ( rr ,j 2 ,...) (10) 

with N' some function of the derivatives of G with respect to L which can be 
deduced from the numerator N . Equation fjlO[) gives a recursion relation for 
calculating the power series expansion of I(k,p). Whenever the absolute value 
of k is not 1, the terms depending on lower order terms of I(k,p) are of the 
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same order than the terms directly proportional to the 7«: the growth of I(k,p) 
is therefore like the one of 7. In the determination of 7, X(k,p) is multiplied 
at least by a 1 , with / the loop order of the correction to the Schwinger-Dyson 
equation, so that such contributions are smaller by a factor 1/n 1 at order n with 
respect to the principal term in 7. 

The scale invariance of quotients G/K with more than 1 edge holds for sums 
of the associated X{ of —2 or less. For a quotient with two loops and three edges, 
there are only three denominators to compensate the two integrations and the 
numerator, which sum to 5, so that the sum of the Xi must be at least —2. The 
quotient with three loops, five edges has at least one unit from the numerator, 
giving again 3x2 + 1 — 5 = 2 to be compensated by the x%. In the case of the 
cube, there is also a four loop quotient, with at least two powers of p 2 in the 
numerator, and eight edges: again the total of the X{ must be —2. The only 
quotients giving terms with k = 1 are therefore the simple edges, giving terms 
proportional to 1(1, 1). Since al(l, 1) is the dominant term in the determination 
of 7 [3J, this will give relative corrections of order 1/n , i.e., of order 1/n 2 
for the three-loop diagram and 1/n 3 for the four loop cube. Corrections of the 
same order in the asymptotic expansion come from the modified order 3 and 
order 4 values of 7 which appear in the recursion for 1(1, 1), viewed as the main 
contributor to 7. 

The situation for the terms with k = — 1 is less simple. However, there is 
certainly one term stemming from the divergence of the full diagram. Indeed, 
there is a pole when the sum of all the x\ is zero, which is compensated when 
taking the derivative with respect to the external momentum, but there are also 
poles for all the negative integer values of this sum, even if the exact numerator 
would be difficult to determine. Now, the number of propagators p is equal to 
3Z — 1 for the term of order I in the Schwinger-Dyson equation. If we now go 
for the recurrence equation for the contribution to 7, a I(—l,p), it becomes 
independent on I, except in the right hand side: 

(-l + 1 (Zad a -l))[a l T(-l,M-l)]=a l N'( 1 , l2 ,...) (11) 

The dominant non-alternating term in the perturbative series for 7, therefore 
receives contributions not only from the order one term in the Schwinger-Dyson 
equation, but also from its corrections at all orders. In any case, the evaluation 
of the coefficient for this asymptotic term should be obtained by a fit to the finite 
number of terms of the perturbative series explicitly computed: the important 
result here is that such a fit should still be based on the same basic series, one 
alternating, the other with alternating signs, and their corrections by powers of 
1/n. 

6 Conclusion 

Our preceding works on the solution to high order of Schwinger-Dyson equa- 
tions were limited to the simplest ones with only one-loop primitive divergences. 
The inclusion of higher- loop primitive divergences seemed dauntless. Indeed, 
the evaluation of the associated Feynman diagrams has to take into account 
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logarithmic corrections of any order in each and every one of the numerous 
propagators in the diagram. Moreover the evaluation of individual contribu- 
tions is but the beginning, since there are numerous possible contributions at a 
given degree in a, corresponding to different ways of obtaining this total degree 
from the sum of the degrees of the individual propagators. Despite all these 
odds, we succeeded to obtain the dominant effect of such a diagram on the 
asymptotic behavior of the perturbative series for the f3— function of the renor- 
malization group. It seems even possible to derive corrections to this dominant 
effect. 

This radical change of perspective stems from the combination of the suc- 
cessive progresses made in our previous works [HOE] with the derivation of 
the general structure of the singularities of the Mellin transform of any diagram 
made in section [3l Here are the successive steps towards this achievement: 

• compute the primitive diagrams with propagators raised to any powers, 
giving a Mellin transform from which the effects of the logarithmic correc- 
tions to the propagators can be obtained through a Taylor expansion [3]. 

• recognize that the asymptotic properties of the Taylor series stem from 
the singularities of the function, so that simple asymptotic formulas for 
these coefficients can be given [2]. 

• for poles of some simple form of the Mellin transform, use the renormaliza- 
tion group to compute the exact contributions stemming from the whole 
Taylor expansion of this pole [3j. For poles with denominators involving 
many variables, this means that we can obtain with a simple operation of 
low computational cost the sum of great number of terms. 

• show that for any diagram, the Mellin transform has only poles of the 
type covered in the preceding item, with computable polynomial residues 
(this work). 

In this article, this whole scheme was illustrated on two specific primitive 
divergences of the Wess-Zumino super symmetric model. We did not completely 
solve the challenge of the evaluation of the numerators, which involves the 
sum of different ones coming from supersymmetry, but we could deduce those 
necessary for the evaluation of the dominant singularities. These partial results 
suggest a high level of symmetry, which should allow a complete description of 
the integral form of the Mellin transform and the computation of all residues: 
we hope to be able to present such results in future works. 

Finally, through the combination of all these insights, we have been able to 
describe the dominant effects of the first corrections to the Schwinger-Dyson 
equations on the asymptotic behavior of the perturbative series. The remark- 
able point is that, even if we cannot compute exactly the full Mellin transforms 
of the primitively divergent graphs, much can be said on the order of magnitude 
of the ensuing corrections without any explicit computation. 

A clear distinction appears between the contributions with constant signs 
and those with alternating signs. The latter ones are dominant in absolute 
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values, have unambiguous sums through the Borel method and receive con- 
tributions which are suppressed by powers of the order from the higher loop 
primitive graphs. It therefore seems possible to have a systematic expansion of 
the higher order contributions, since only a finite number of primitive graphs 
give contributions to each order of the expansion. This could be trumped by the 
number of primitive divergences, which is a finite fraction of the total number 
of graphs, but the result should in any case hold in the large N limit where only 
planar graphs contribute. If we had only these terms, we could think that a con- 
structive determination of the renormalization function is at hands. However 
the parts with constant sign are there and cumulate difficulties: they are asso- 
ciated to singularities of the Borel transform for positive arguments, implying 
ambiguities in the summation process, and every order of the asymptotic ex- 
pansion receives contributions from all possible terms in the Schwinger-Dyson 
equation. 
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